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ABSTRACT 



The dissipation method is used to obtain estimates for 
the friction velocity, U*, as well as values for the neutral 
drag coefficient, for data collected from a coastal 
tower off San Diego, California. Cp^j is found to be indepen- 
dent of the ten-meter height windspeed, Uj^q, for velocities 
between 4-9 m/sec. Its value is estimated to be 
(0.94+0.4)10^ which compares well with values by Smith 
(1980) and Large and Pond (1981). Definite trends in Cjjj| 
with fetch and sea state are also observed. Drag coeffcient 
estimates are found to be higher for short fetch than for 
long fetch conditions. Cjjjq- is also seen to increase sharply 
just before frontal passages and during sea breeze condi- 
tions when the waves are actively growing. With the wind- 
speed and wave field reaching equilibrium, ^DN is found to 
decrease with time to a smaller and more constant value. 
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I. INTRODUCTION 



The Marine Atmospheric Boundary Layer (MABL) is the 
region adjacent to the air-sea interface where the vertical 
variation of the mean stresses is large compared to the 
horizontal gradients (Holton, 1979). In this layer, turbu- 
lent shear stresses due to both molecular properties 
(viscosity and diffusion) and turbulent eddies coexist. 
These eddies are very effective mixing agents which serve to 
transfer momentum, heat and other quantities to and from the 
earth's surface, at a rate much faster than that for molec- 
ular diffusion. Frictional forces due to molecular viscosity 
are usually neglected by means of scale analysis (Holton, 
1979), and the emphasis has been placed on parameterizing 
the fluxes resulting from the turbulent eddies. 

The mean vertical fluxes of momentum, sensible heat and 
latent heat are (Fleagle and Businger, 1980): 



where t, H and E are the momentum flux (wind stress), 
sensible heat flux and latent heat flux, respectively. u’ 
and w' are the horizontal and vertical turbulent wind compo- 
nents, while 0 ' and q' represent the turbulent fluctuations 



latent heat of vaporization, is the air density and Cp is 
the specific heat of air at constant pressure. 

In order to estimate these fluxes using routine surface 
layer observations, bulk transfer formulae were developed 
(see Fleagle and Businger, 1980): 



T 

H 



E 




(la) 

(lb) 

(l c) 



of potential temperature and specific humidity. L£ is the 



T 




(2a) 
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H 



E 



- UJ(T^ - T^) 
PaLfCE(Uz - U^)(qo - ^z) 



(2b) 

(2c) 



where and represent the windspeed, temperature and 

humidity specified, by convention, for a reference height, 
z. Uq, Tq and q^ are the same quantities measured at the 
surface. Cp, and Cj- are dimensionless quantities repre- 
senting bulk transfer coefficients for momentum, sensible 
heat and latent heat, and are for the height, z. These are 
known as the drag coefficient, the Stanton Number, and the 
Dalton Number, respectively. 

The turbulent exchange processes greatly influence the 
general circulation of the atmosphere and thus have tremen- 
dous impact on countless social and environmental concerns . 
Central to their understanding is a determination of the 
nature of the bulk transfer coefficients, particularly the 
drag coefficient, C^. This quantity encompasses the mecha- 
nism responsible for the wind stress at the surface. Since 
most of the parameters used to describe the bulk formulae 
are routinely measured, ascertaining a reliable value for 
based on these quantities is essential for accurately 
defining the surface stress, t- Further, Cp is an important 
factor driving many of the current ocean wave, acoustic and 
boundary layer models (Geernaert, 1985a). 

Considerable effort has gone into examining the over- 
water drag coefficient. It has been shown to be windspeed 
and stability dependent, and may also vary with factors such 
as fetch and depth (Geernaert, 1983). Studies also indicate 
that the overwater drag coefficient tends to have a magni- 
tude on the order of 10 . Its value exhibits considerable 

scatter under varying environmental conditions. Higher 
values are usually associated with an unstable atmosphere, 
short fetch and/or shallow seas, implying a wave dependence 
(Geernaert, 1983). The more recent use of remote sensing to 
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study winds over the ocean requires an understanding of this 
wave dependence since it assumes that the wave structure is 
related to the wind stress. Clearly, the processes 
affecting the drag coefficient merit considerable further 
study. This thesis will investigate the dependence of the 
drag coefficient on windspeed, fetch and wave state. 

An important scaling parameter, based on the wind 
stress, is the friction velocity , U*, which may be defined 
as : 

T = (3) 

The quantity U,v is related to most air- sea exchange 
processes such as heat and moisture fluxes (Businger, 1973), 
the surface drift (Hicks, 1972) and ocean mixed layer 
dynamics (Kraus, 1972). 

Data considered in this investigation were from the 
Tower Ocean Wave and Radar Dependence (TOWARD) experiment. 
Located offshore of Mission Beach at San Diego, California, 
the Naval Ocean Systems Center (NOSC) tower provided a 
stable platform with a capability for mounting a variety of 
sensors to measure the variables needed to parameterize the 
bulk formulae. Naval Postgraduate School (NFS) personnel 
made atmospheric surface layer measurements of both mean and 
fluctuating velocity, temperature and humidity. The NPS 
measurements are listed in Table I. Wave information based 
on radar signal returns was collected through investigations 
from the Naval Research Laboratory, the Jet Propulsion 
Laboratory and the University of Kansas. The Scripps 
Institute of Oceanography measured tower frequency wave 
spectra with an array of sub-surface pressure transducers. 

The measurement strategy for the TOWARD Experiment was 
originally composed of two intensive data collection 
periods; September/October 1984 and January/ February 1985. 
These periods were chosen because of their differing average 
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windspeeds and resultant sea states. With its relatively low 
average windspeeds, the fall measurement period had calm to 
moderate sea state conditions. The winter period, with its 
higher average windspeeds, was designed to focus on more 
active sea states. Equipment problems, however, forced the 
postponement of Phase I to October/November 1984 and Phase 
II to March/April 1985. 

Figs, 1.1 and 1.2 provide two views of the NOSC tower. 
Located in 17-18 meters water depth, the tower provided a 
stable platform for data collection. A structural extension 
of the tower from the south side allowed wave rider sensors 
to be placed away from the main tower structure. 

Measurements of wind speed and direction, wave height, 
tide level and near surface currents were recorded routinely 
during intensive measurement periods. Air and sea-surface 
temperatures were also recorded, along with wind stress data 
that was obtained with a sonic anemometer and a hot film 
sensor . 
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TABLE I 




NFS Measurements 




Measurement 


Sensor 


Sea-Surface Temperature 


Floating platinum 
thermometer 


Mean Surface Layer: 




Wind (speed and direction) 

Temperature 

Humidity 

Radiation (solar and IR) 


Cup anemometer, 
vane 

Aspirated platinum 
thermometer 
Dew cell 
(cooled mirror) 
Eppleys (up and 
down) , Total 


Turbulent Surface Layer: 




Wind stress: direct 
Wind stress: indirect 
K. E. dissipation rates 
Sensible heat flux 

Latent heat flux 

Spectra: windspeed 
temperature 
humidity 


Sonic anemometer 
Hot films 
Hot films 

Platinum thermometer 
and sonic anemometer 
Lyman- alpha and 
sonic anemometer 
sonic anemometer 
Platinum thermometer 
Lyman- alpha 
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Figure 1.1 The NOSC Tower. 
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Figure 1.2 The NOSC Tower. 
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II. REVIEW OF EXISTING SURFACE LAYER THEORY 



The surface layer comprises approximately the lowest ten 
per cent of the MABL and is defined as the layer of air 
nearest the ground where atmospheric parameters may be 
scaled with surface fluxes. Within this layer, an approxima- 
tion adopted is that the fluxes are independent of height 
(Businger, 1973). 

In Eqn. (2a), the surface speed, , is a parameter 
necessary for estimation of the drag coefficient. However, 
because is very small compared to U^q, Uq is assumed to 

be zero in the bulk definitions. With this assumption, Eqns . 
(2a) and (3) may be combined into the form: 

Cdio = (U*/Uio)^ (^) 

Note that the drag coefficient in Eqn. (4) is defined for a 
height of 10 meters and is denoted as 

The similarity theory which described the mean surface 
layer wind profile is (Monin and Obukhov, 1954): 

d^ldz = (U,v/kz)(^^(z/L) (5) 

where 5U/^z represents the vertical gradient of the wind 
speed and k is the von Karman constant, which ranges in 

value from .35 to .43 (Panofsky and Dutton, 1984). This 
study assumes a value of 0.4 for this constant. is 

a dimensionless stability function. One parameterization of 
this function is that presented by Businger, et . al., 
(1971): 

= [1 - a( 2 /L)]-l/^ (z/L<0) (6a) 

<^ni = [1 ^ j3(z/L)] (z/L>0) (6b) 
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where a and j3 are diabatic profile constants, with approxi- 
mate values of 15±.l and 5i0 . 5 , respectively (Geernaert, 
1983). z is the measurement height and L is the 
Monin-Obukhov length, defined as the height at which the 
buoyant energy production equals the shear production 
(Businger, 1973). z/L is defined as: 

z/L = -gkzw'T^'/TQ’ (U*3) (7) 

The virtual temperature fluctuation, T.^.' , is related to 
the air temperature fluctuation, T', and specific humidity 
fluctuation, q', by: 

T^' = T’ (1 + .61q) + . 61Tq' (8) 

Thus : 

w’T^' = w^'(l + .61q) + .61Tw^’ (9) 

where g is the gravitational acceleration. The first term 
on the right hand side of Eqn. (9) represents the tempera- 
ture flux and the second term describes the humidity flux. 
These terms are further estimated by bulk equations where: 

~ = ChUio(To - T^) (10a) 

“'q' = - q^) (10b) 

Both Cg and C^j are on the order of 10 , but are less well 

defined than the drag coefficient (Geernaert, 1983). 

Under neutral or adiabatic conditions, 

■#>m(a/L) = 1 

and Eqn. (5) reduces to: 

3U/<9z = U-v/kz 

Integrating Eqn. (12): 

.o(z) * 

/u(ro') = (U*/k)4^dz/z 

yields the well known logarithmic wind profile: 



( 11 ) 

( 12 ) 
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u 



(U*/k)ln(z/zQ) 



(14) 



where the roughness length, z^ , is a measure of the effec- 
tive surface roughness. 

For non-neutral (diabatic) conditions, the stability 
parameter must be included. Integrating Eqn. (5) yields the 
stability dependent logarithmic wind profile (Paulson, 
1970) : 

U = (U*/k)[ln(z/zQ) - ,//] (15) 

where \(/ represents the stability parameter. 

For unstable stratifications (z/L<0): 

yp = 21n[(l+x)/2] + ln[(l+x^)/2] - 2(tan)"lx + (n/2) (16) 

where : 

X = [(1 - a(z/L)]^/'^ (16a) 

For stable conditions (z/L>0): 

xp = - j3(z/L) (17) 

The logarithmic wind profile is important primarily for 
two reasons. First, it can be used to estimate the ten 
meter height wind speed, U^q. Beginning with Eqn. (14), the 
wind speed at a particular height, z, may be described as: 

= (U-.v/k)[ln(z/z^) - i//^] (18) 

At ten meters elevation: 

Uio = (u*)/k[ln(10/z^) - (19) 

Combining Eqns . (18) and (19) yields the following relation- 

ship : 

Uio = > (u*/k)[ln(10/z) - (20) 

This particular equation is important because it is not 
always possible to place sensors on ships or platforms at 
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can be 



the required ten meter height. As a result, 
estimated regardless of the height of the instrument. 

Second, the log wind profile is important in determining 
the neutral drag coefficient, which may be defined as 

the wind stress exerted on the surface for neutral stratifi- 
cations. Estimation of the neutral drag coefficient is 
important because it eliminates variations in Cjj parameteri- 
zations due to differing stability conditions. Many 
modellers prefer to be provided with instead of 

because the in-situ bulk weather data may be used to esti- 
mate a stability correction. Eqn. (24) may then be applied 
to a given Cj^jj and \p . 

A. GENERAL DRAG COEFFICIENT REPRESENTATIONS 

Eqns . (4) and (15) yield: 

Cd = (U*/U2)^[k/[ln(z/Zo) - i//^(z/L)]]2 (21) 

Under neutral conditions, where \p=0: 

= [k/ln(z/zQj^)]^ (22) 

where Cj^j^ and z^^^ are the drag coefficient and roughness 
length under neutral conditions. Rearranging Eqn. (22) 
yields : 

= lri(z/zQj^) (23) 

Combining Eqns. (21) and (23) produces the following general 
result (Geernaert, 1983): 

Cd = * (l/k)ln(z^^/zQ)]-2 (24) 

The drag coefficient is known to vary under differing 
environmental conditions. Studies have indicated that C^^ 
increases more rapidly with increasing winds over water than 
over land (Hsu, 1974). Geerneart (1983) discussed the 
application of Eqn. (24) to both land and water cases. Over 
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land, the roughness length is assumed to be largely indepen- 
dent of the surface stress. Thus, = z^ and Eqn. (24) is 
reduced to: 

Cp = (^DN ~ ^ (25) 

Over water, both the drag coefficient and, hence, the 
roughness length show considerable windspeed and stress 
dependence (Kitaigorodskii , 1973). Charnock (1955) was 
among the first to recognize that Cjj increases considerably 
with windspeed over the sea. He postulated that the rough- 
ness length, and hence, the neutral drag coefficient was 
proportional to the wind stress through the following 
scaling argument: 

Zo = a'U*^/g (26) 

where «' is the Charnock constant with a typical value of 
.015-. 020 (Geernaert, 1983). Charnock' s work served as a 
basis for subsequent parameterizations of the drag coeffi- 
cient and the roughness length in terms of sea state 
elements . 

Kitaigorodskii (1973) further described z^ based on wave 
spectral densities. Since that time, z^ and Cp have been 
shown to be dependent on a variety of other factors that are 
linked to the character of surface waves such as fetch, 
depth and wave age. Cp , therefore, can be evaluated as: 

Cp = *^dN^^ 10’ waves) (27) 

As a result, a different relationship is often used to esti- 
mate the drag coefficient over the ocean. • Applying 
Charnock' s relation to Eqn. (24) yields (Geernaert, 1983): 

% = ^ (l/k)ln(Cj5N/CD)]‘2 (28) 
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B. DETERMINING THE WIND STRESS BY THE DISSIPATION METHOD 



Since estimation of the drag coefficient requires values 
for the friction velocity, it is important to examine how 
the momentum flux is measured. Two primary methods for 
obtaining t are the eddy correlation technique and the 
dissipation method. The former method involves measurement 
of the turbulent wind components u', v’ and w' , and ideally 
requires a steady platform. Because of the low frequency 
motions inherent in sea-going vessels on which most marine 
studies are conducted, the dissipation method is often 
substituted for oceanic wind stress experiments. 

The dissipation method is an indirect approach in which 
the rate of dissipation of turbulent kinetic energy is used 
to determine U*. This method requires application of the 
turbulent kinetic equation: 

d'^/dt = U*^5U/(9z+ (qw'T^' )/T q-5/(9z[w' e+ (l/p^)w'p' ] -G (29) 

where d^/dt is the local time rate of change of turbulent 
kinetic energy (TKE), Uvf^5U/(9z is the wind shear produc- 
tion, (qw'T^')/!^ represents the buoyancy forces, 
3/6z[w'e + (l/pa,)w'p'] is the redistribution of the flux of 
TKE and work done by pressure fluctuations and g is epsilon, 
the rate of dissipation of turbulent kinetic energy. 

Steady state conditions are assumed. The divergence term 
is small relative to the remaining terms in the equation and 
is therefore neglected (Large and Pond, 1981). Eqn. (29) 
then simplifies to: 

u-^au/az + [(g^r^')/T^] - g = o ( 30 ) 

After combining Eqns . (5), (7) and (30), and rear- 
ranging, an expression for the friction velocity is 
obtained : 

U* = [Gkz/ [^jjj(z/L) - (z/L)]]^/^ (31) 
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For unstable conditions, the denominator is often replaced 
by : 



<^^(z/L) = <^>ni(z/L) - (z/L) (32) 

Eqn. (31) may then be written as: 

U* = [£kz/,^^(z/L)]l/3 (33) 

with stability corrections for dissipation scaling (Wyngaard 
and Cote ' , 1971 ) : 

= (1 + 0.51z/L|^/^)^/^ z/L<0 (34a) 

= (1 + 2.5|z/L|^/3)^/^ z/L> 0 (34b) 



The friction velocity, U*, can therefore be estimated in 
terms of the dissipation rate, £. Kolmogoroff theory is then 
employed to determine a value for e based on high frequency 
wind fluctuations. All spectra can be divided into three 
portions: the energy producing subrange, the inertial 
subrange and the dissipation subrange (Panofsky and Dutton, 
1984). Large eddies from the energy containing subrange 
cascade kinetic energy and momentum down to smaller eddies, 
eventually reaching the dissipation region where molecular 
viscosity converts the kinetic energy to heat. Between these 
two regions is the inertial subrange. In this region, energy 
is neither produced nor dissipated. Kolmogoroff theory 
states that in the inertial subrange of isotropic turbu- 
lence, the high frequency fluctuations of the wind velocity, 
u' can be expressed by a one-dimensional power spectrum 
(Schacher, et . al., 1981): 

Su(k*) = (35) 

where Su is the energy density spectrum of the horizontal 
windspeed, is a coefficient of magnitude 0.50 (Champagne, 
et. al., 1977), and k* is the wavenumber. 
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Applying Taylor's frozen turbulence hypothesis: 



k* = 2irf/U (36) 

equation (35) becomes: 

Su(k*) = a*€^/^(2irf/U)-^/3 (37) 



The Kolmogoroff Spectrum can now be integrated between 
two wavenumbers, designated and k.,vj^ to give the energy, 

^^Ak-'- iriterval (Schacher, et . al., 1981): 



<^^Ak-. = Su(k*)dk* 

= 1. 5avf€^^^ 



Combining Eqns • (36) and (38) yields: 

= 0.230(£0)2/3(f^-2/3 _ fj^-2/3) 






(38) 

(38a) 



(39) 



where f (k^v^^U) /2 and f j^= (k,vj^U) /2 . Rearranging Eqn. (39) by 
using Eqn. (35) produces: 



U.,v = 2.65auAf [U<^g(z/L)]]^/^(f 3 ^"^/^ - 



(40) 



For the TOWARD Experiment, = 5 Hz and = 50 Hz, and 

z=22 meters. Substituting these values reduces Eqn. (40) to: 



U* 



7.81au^f[U<^g(z/D] 



-1/3 



(41) 



Because 0^(z/L) also depends on U*, Eqn. (41) must be 
solved iteratively. 



C. IN- SITU HOT FILM CALIBRATION 

Calibration of the hot film sensor was accomplished by 
recording the changes in voltage produced by the winds 
passing over the sensor. The relationship may be expressed 
by: 

+ BUrel^/^ 
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where V is the voltage, is the plotted calibration inter- 
cept and Urel is the relative windspeed. B represents a 
calibration factor that must be continuously determined for 
each data collection period. 

B may be estimated either through the dynamic or the 
regression methods. The dynamic method begins by differen- 
tiating Eqn. (42): 

B = 4VUrel^/^(dV/dU) (43) 

Over each data collection period, dV/dU is estimated by 
av/c7u> where ay and ay are the standard deviations of U and 
V, respectively. Averaging V and Orel and solving for B 
yields : 

B = 4VUrell/^(ay/au) (44) 

Rearranging Eqn. (44) and substituting into Eqn. (40) 
produces the following estimate for U*: 

U* 

2.81[(4V/B)Urell/2]jy^^^[2/U(^^(z/L)]l/3(fl'2/3 - ( 45 ) 

The regression method is based on values of V and Urel 
over several time periods with varying average relative 
windspeeds . V is then usually plotted as a function of Urel. 
According to Eqn. (42), the points should fall along a line 
with a slope, B, where represents the y-intercept. 

Measurement-relative scatter in both Urel and V can 
cause errors in the slope. Therefore, it is sometimes advan- 
tageous to choose V as the independent variable if it shows 
less scatter. In cases where it is difficult to determine 
which variable is the least reliable, two regressions are 
calculated, with Urel and V alternated as independent vari- 
ables. B is then determined from the average of the two 
slopes . 
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D, WIND-WAVE COUPLING DYNAMICS 



The total wave field energy is divided into swell and 
wind wave components: 

= a^s ^ ty^ww (^6) 

n 

where a i represents the total energy in the wave field 
2 2 

and a ^ and cr are the energies contained by the swell and 
wind waves, respectively, a cannot be predicted based on 
local scaling, o has received considerable study and a 
great deal of theory exists concerning this parameter. The 
approach currently is to describe in terms of a local 

equilibrium spectrum. This theory assumes that wind forcing 
controls the dynamics of every wave component of the wind 
wave spectrum. Each wave is characterized by an angular 
frequency, n', where: 

n' = lirt (47) 

Gravity is the primary restoring force for these waves . 

The local equilibrium energy density spectrum, S(n'), is 
defined as : 

= /rs(n')dn’ (48) 

and S(n') may be parameterized in terms of Phillips scaling 

arguments (Phillips, 1980): 

S(n') = j3g(n’)-^ n’>nQ’ (49) 

Beta is the Phillips coefficient which can vary from .006 to 
.020 (Geernaert, 1983). n^ ' represents both the wind-wave 
peak frequency and the boundary between the swell and equi- 
librium ranges. Waves of frequency n'<nQ' represent swell. 
Wind waves are defined to be those with frequencies of 

n'>nQ'. n^ ' is usually approximated by: 

%' = B/Uio (50) 
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for steady state conditions (Kraus, 1972). 

Beta is often described in terms of dimensionless fetch, 

A# 

x: 

X = gx/(u*^) (51) 

where x is the upwind fetch. Large values of jS are associ- 
ated with small values of fetch. Beta has also been parame- 
terized over an applicable range of .25 to 1.0 hz by 
(Geernaert, et.al., 1985a): 

j3 = 0.005 + .002n' + 1.5 (Cq/U*)'^ (52) 

where represents the phase speed of the largest wind 

wave. (CQ/Uv^) is defined as the wave saturation parameter, 
and this parameter behaves in the same manner as the wave 
age C/U. (Cq/U-v) is a dimensionless indicator of the degree 
to which the wind waves have reached steady state and is 



frequently related to the non-dimensional fetch as (Wu, 
1985): 

(Cq/U-/) = 0.05[gx/(U,v^)]3/10 (53) 

Hsu (1974) related the wave saturation parameter to the 
roughness length, , by defining the Charnock coefficient 
in terms of the slope of the wave of frequency n^ ' . His 
formulation became: 

Zq = (H/2^)(Cq/U*)'^ (54) 



where H is the significant wave height . z^ can be seen to 
vary according to the local wind wave spectral peak. 

Numerous wave models have incorporated these relation- 
ships in an attempt to relate variations in the surface wave 
energy spectrum to variations in the magnitude of the drag 
coefficient. Davidson (1974) formulated an empirical model 
in which the drag coefficient is represented by: 

Cj3 = [k/[ln(z/Zo) - ^ ^ b[(Co/U*) - 26.3]]]2 (55) 
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where b is an empirical constant found to be 0.13. 

Others have employed indirect methods to calculate Cj^ 
using roughness length parameterizations . Byrne (1982) 
derived a roughness length in the form of an integrated 
slope spectrum: 

2o = CB/J^S(n’)(n')2dn' (56) 

where Cg is an adjustable coefficient. 

Kitaigorodskii (1973) derived a roughness length model 
in the following form: 

Zq = 2/30(/J^S(n')e-2kg/n’U*dn’)l/2 (57) 

One purpose of the TOWARD Experiment was to evaluate 
scatterometer return signals in order to examine surface 
roughness as a predictor of windspeed over the ocean. These 
signals can be correlated with windspeed since they increase 
with increasing surface roughness. The surface roughness is 
dependent on t, which from Eqn. (2a) is a function of and 
the ten meter height windspeed, scatterometer, 

however, evaluates U^q based on signal reflection from the 
surface (Geernaert , 1983), while the drag coefficient has 

been treated as a constant or as a simple function of U^q. 
Thus, the dependence of the drag coefficient on the surface 
roughness is not fully understood. Because wave growth and 

decay affect the magnitude of (Geernaert, 1985a), scat- 

terometer predictions of windspeed based on the bulk 
formulae could be improved by including roughness as an 
additional variable. 

The TOWARD Experiment was designed, in part, to look for 
a trend in Cj^ as a function of time over the course of the 
sea breeze. The sea breeze is characterized by growing winds 
and sea state along with corresponding increases in surface 
roughness. By observing the behavior of this coefficient 
with time, the interaction between wind and wave state may 
be better understood. 



29 



III. SYNOPTIC METEOROLOGY DURING THE TOWARD EXPERIMENT 



The TOWARD experiment was conducted at the Naval Ocean 
Systems Center (NOSC) tower which is located two kilometers 
offshore of Mission Beach, San Diego. Minimal surface layer 
climatological data is available for the site. Blanc (1981) 
summarized the environmental conditions for San Nicholas 
Island (SNI), California, which is located roughly 150 km 
northwest of the NOSC tower. While wind and wave conditions 
are more vigorous at SNI than those at the tower, the air 
and sea-surface temperatures are consistent between the two 
sites. Blanc determined the most frequent true wind direc- 
tion to be from the northwest throughout the year with mean 
windspeeds of 7.7 m/s in March and 6.2 m/s.for both October 
and November. During the months of March, October and 
November, when the TOWARD experiment was conducted, the mean 
climatological daily temperatures (°C) are 15.8 for March, 
20.4 for October and 18.4 for November, with average daily 
minimum temperatures of 12.4, 16.0 and 14.9, respectively. 
Mean climatological sea-surface temperatures for the same 
months are 14.4 for March, with values of 18.2 and 16.6 for 
October and November. 

The mean conditions for San Nicholas Island differ some- 
what from those reported onshore at Lindbergh Field in San 
Diego. The mean daily maximum temperatures for March, 
October and November are 18.9, 23.3 and 21.2, respectively, 
while the average daily minimum temperatures are 10.1, 14.7 
and 10.8. In March and October, the prevailing winds are 
generally west-northwesterly but shift to a northeasterly 
direction in November. Mean windspeeds are also smaller 
with values of 3.3, 2.9 and 2.5 m/sec reported for March, 
October and November, respectively. 
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The TOWARD Experiment was designed to cover a large 
range of wind conditions, particularly for high winds. 
Measurement periods were therefore chosen to be from 
mid-September through mid-November 1984 (Phase I) and March 
1985 (Phase II). 

The almost daily presence of the local sea breeze is a 
factor in the environmental conditions at the NOSC tower. 
Sea breeze circulations develop when air over land is heated 
much more than the air over the ocean. Since the land 
usually becomes warmer than the sea during the day, density 
differences over the two surfaces bring cool air from the 
sea toward the warmer land at low levels . A return flow of 
warm air from the land to the ocean occurs at upper levels. 
At night, radiative cooling takes place over the land while 
the sea surface maintains a more uniform temperature due to 
the large heat capacity of the ocean. The result is 
nocturnal, low level flow from the land toward the sea 
(Hess, 1979). 

Because the largest horizontal density gradient occurs 
at the coastline, the sea breeze dominates in the coastal 
zone, while the geostrophic wind (V„) becomes a more impor- 
tant factor for describing surface layer energetics farther 
offshore. It follows that waves generated by the sea breeze 
would tend to be of relatively short fetch. Fig. 3.11 from 
31 October provides an example of synoptic conditions for 
Southern California during the TOWARD Experiment. The large 
scale flow (Vg) can be seen to originate from the northwest. 
Ship observations, 300-400 km offshore, also indicate a 
northwesterly flow. These wind observations from offshore 
sites show little diurnal variation. At 1000 local time, 
however, reports from Lindbergh Field, San Diego, show winds 
from a south-southwesterly direction, indicating smaller 
scale, local- forcing nearer the coast. A significant diurnal 
shift in wind direction can also be seen in these reports. 



31 



Thus, the sea breeze appears to exert an influence less than 
300 to 400 km from the shore. 

Since the sea breeze is of a different scale than the 
synoptic, geostrophic wind, computation of Vg provides an 
indication of the relative magnitudes of these two types of 
forcing near the NOSC tower. Calculations of the geos- 
trophic wind from synoptic charts were made for 1000 local 
(1800 GMT) and 1600 local (2400 GMT). 

Vg = (l/fp^)(5p/an) (58) 

where Vg represents the geostrophic wind, f is the coriolis 
parameter, p is the atmospheric pressure and n is a unit 
vector normal to the isobars. These times were chosen in 
order to compare the differences in sea breeze magnitudes 
between morning and afternoon. Typically, the strongest sea 
breeze occurs in mid-afternoon when the land surface has 
heated to its daily maximum. Since final surface charts 
were not available for 2100 GMT (1300 local), analyses for 
2400 GMT were used instead. These values were compared with 
windspeed observations collected at the tower and with data 
from ship reports in the Southern California area. 

Data collected at the tower indicated a general increase 
in windspeed from 1000 to 1600 local (1800 GMT to 2400 GMT). 
The geostrophic wind velocities based on Eqn. (58) and the 
synoptic pressure maps were smaller and showed much less 
diurnal variation than those of the tower (V^.). Ship obser- 
vations from 300-400 km offshore indicated consistently 
higher windspeeds for both time periods, usually ranging 
from 5-15 kt . Although the computed V and ship reported 

O 

surface winds (Vg) showed considerable differences in magni- 
tude, the wind directions from offshore sites during the 
experiment showed no statistically significant diurnal 
effects. Tables II and III show these results. 



32 













TABLE II 








Computed Vg & 


Ship 


Reported 


Sfc 


Winds 


(Phase 


I) 








( Vf , V 


g Vg 


in m/ s ) 






Date 


Time 


Vt 


Dir 




Dir. 


Vs 


Dir. 


D. (km) 


10/2 


18 


Z 


2.0 


87 


<1 


NW 


2.6 


N 


400 




24 


Z 


— 


— 


1.2 


NW 


7.7 


N 


400 


10/3 


18 


Z 


2.5 


297 


1.4 


NW 










24 


Z 


5.3 


309 


<1 


NW 


y.i 


N 


400 


10/6 


18 


Z 


— 


— 


<1 


NNW 


2.6 


W 


400 




24 


Z 


5.9 


313 


1.4 


NNW 


7.7 


N 


400 


10/7 


18 


Z 


4.6 


335 


<1 


NW 


2.6 


N 


400 




24 


Z 


- 


— 


1.7 


NW 


5.2 


N 


400 


10/11 


18 


Z 


— 


— 


1.4 


NW 


7 . 7 


N 


400 




24 


Z 


— 


— 


1.1 


NW 


7.7 


NW 


400 


10/12 


18 


Z 


1.0 


154 


1.0 


NE 


10.2 


N 


400 




24 


Z 


8.3 


317 


1.4 


NE 


10.2 


N 


400 


10/15 


18 


Z 


6.0 


180 


1.4 


NE 


10.2 


N 


400 




24 


Z 


5.4 


293 


1.2 


NE 


12.8 


N 


400 


10/16 


18 


Z 


— 


— 


<1 


NNW 


2.6 


NW 


400 




24 


Z 


5.4 


203 


<1 


NW 


7.7 


N 


400 


10/17 


18 


Z 


6.7 


280 


1.4 


N 


10.2 


N 


400 




24 


Z 


6.2 


181 


1.7 


N 


7.7 


NW 


350 


10/18 


18 


Z 


— 


— 


1.0 


ENE 


10.2 


NNW 


400 




24 


Z 


6.8 


337 


<1 


ENE 


7 . 7 


N 


400 


10/20 


18 


Z 


5.6 


255 


<1 


N 











24 


z 


6.0 


267 


1.4 


NW 


10.2 


N 


400 


10/21 


18 


z 


1.8 


127 


<1 


NNE 


— 








24 


z 


— 


— 




NW 


7 . 7 


N 


400 


10/23 


18 


z 


— 


— 


<1 


NE 


5.1 


N 


400 




24 


z 


8.7 


288 


<1 


ENE 


— 


— 


— 


10/25 


18 


z 


— 


— 


<1 


SE 


2.6 


N 


400 




24 


z 




— 


<1 


SE 


2.6 


NW 


400 


10/27 


18 


z 


— 


— 


<1 


N 


— 




— 




24 


z 


4.4 


221 


1.7 


N 


5.2 


N 


400 


10/28 


18 


z 


— 


— 


<1 


NE 


— 




— 




24 


z 


3.9 


298 


<1 


NW 


7.7 


N 


400 


10/29 


18 


z 


— 


— 


1.1 


NE 










24 


z 


4.4 


315 


<1 


NE 


i ’2 


NW 


400 


10/30 


18 


z 


1.4 


231 


<1 


NNW 


— 


— 


— 




24 


z 


3.7 


294 


2.8 


NW 


7.7 


N 


400 


10/31 


18 


z 


2.8 


218 


1.7 


NW 


— 


— 


— 




24 


z 


5.3 


296 


1.7 


NW 


— 


— 


— 


11/2 


18 


z 


2.6 


298 


1.7 


NNW 


7 . 7 


N 


400 




24 


z 


6.8 


325 


1.4 


NW 


7 . 7 


N 


400 


11/3 


18 


z 


— 


— 


<1 


NW 


2.7 


SW 


400 




24 


z 


3.9 


289 


1.4 


NE 


10.2 


N 


400 


11/4 


18 


z 


1.1 


170 


<1 


SE 


7.7 


N 


400 




24 


z 


2.8 


306 


<1 


SE 


— 


— 


— 


11/5 


18 


z 


1.4 


158 


<1 


N 


5.1 


N 


400 




24 


z 


1.5 


199 


2.1 


NW 


7 . 7 


N 


400 


11/7 


18 


z 




— 


< 1 , 


NW 


7 . 7 


N 


400 




24 


z 


4.1 


257 


1.4 


NW 


5 . 1 


N 


400 
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TABLE III 



Computed Vg & Ship Reported Sfc Winds (Phase II) 









(V 


t’ Vg 


and Vg 


in m 


/s) 






Date 


Time 


Vt 


Dir . 


Vg 


Dir. 


Vs 


Dir . 


D. (km) 


3/12 


18 


Z 


_ _ « 


... 


1.4 


NNW 


5.1 


N 


400 


24 


Z 


6.0 


260 


<1 


NW 


10.8 


N 


400 


3/13 


18 


Z 


— 


— 


<1 


NE 


7.7 


NW 


250 














5.1 


NW 


275 




24 


Z 


3.0 


260 


1.9 


N 


7.7 


N 


400 


3/15 


18 


Z 


— 


— 


<1 


NW 


5.1 


N 


400 


24 


Z 


4.8 


415 


<1 


NW 


7.7 


N 


400 


3/16 


18 


Z 


4.0 


003 


<1 


NW 


2.6 


NW 


400 


24 


Z 


4.9 


340 


<1 


NW 


7.7 


N 


400 


3/18 


18 


Z 


— 


— 


1.4 


NW 


7.7 


N 


400 


24 


Z 


9.4 


296 


2.4 


NW 


10.3 


N 


400 


3/19 


18 


Z 


— 


— 


2.8 


NW 


10.3 


N 


400 


24 


Z 


4.8 


275 


1.7 


NW 


12.9 


N 


400 


3/20 


18 


Z 


— 


— 


1.3 


NNW 


7.7 


NW 


400 


24 


Z 


5.9 


312 


1.4 


NW 


7 . 7 


N 


400 


3/21 


18 


Z 


— 


— 


2.1 


NW 


10.3 


N 


400 


24 


Z 


3.2 


315 


1.8 


NW 


18.0 


N 


400 


3/22 


18 


z 




— 


1.4 


N 


12.9 


N 


400 


24 


z 


5’6 


322 


1 


NE 


12.9 


N 


400 


3/24 


18 


z 


— 


— 


<1 


NE 


5.1 


NW 


275 














5.1 


N 


400 




24 


z 


— 


— 


<1 


NNW 


7.7 


N 


400 


3/25 


18 


z 


— 


— 


1.4 


NNW 


7 . 7 


N 


400 


24 


z 


6 . 1 


220 


1.7 


NNW 


10.3 


N 


400 


3/26 


18 


z 


3.6 


218 


1.7 


NW 


7 . 7 


N 


400 


24 


z 


3.6 


265 


<1 


NW 


10.3 


N 


400 


3/27 


18 


z 


7 . 6 


233 


<1 


NW 


5.1 


NNW 


400 


24 


z 


5 . 5 


264 


2.1 


NW 


7.7 


NNW 


400 


3/28 


18 


z 


7.5 


235 


2.8 


NW 


12.9 


N 


400 


24 


z 


— 


— 


2.8 


NW 


10.3 


N 


400 



A. THE TOWARD SYNOPTIC AND SURFACE LAYER CONDITIONS 

This section will be a presentation of the surface layer 
and larger scale meteorological conditions for the TOWARD 
experiment. The data set is divided into two sections, Phase 
I and Phase II. Phase I covers the period from 2 October to 
8 November 1984, while Phase II covers the period of 12 to 
28 March 1985. Each phase is discussed separately. 
Individual surface analyses for 1800 GMT will be presented 



34 



in the text for days on which data were collected. 
Additional figures representing wind speed and direction, as 
well as air and sea- surface temperatures can be found in 
Appendix A. The latter figures present time series of 
discussed quantities and provide a reference for the 
following discussion. Unless otherwise noted, the times are 
local (PDT for Phase I and PST for Phase II). 

1. Phase I 

During the first week of the TOWARD experiment, 
winds were from the northwest and windspeeds were light to 
moderate. A maximum value of roughly 7 m/s (Fig. A.l) was 
recorded on October 6th. Strong local diurnal variations due 
to the sea breeze also prevailed in the first week, which is 
particularly evident in the data from 3 and 6 October. On 
those dates, windspeeds increased Significantly from 0900 
PDT, remained steady for several hours, and then decreased 
markedly throughout the sampling period which ended at 2215. 

Data collection resumed at 1235 PDT on October 6th. 
Recorded windspeeds remained fairly steady from this time to 
1500, when values began decreasing. Data were only collected 
for one hour on 7 October. Windspeeds increased steadily 
from about 4 to 6 m/s over this period. 

Synoptic charts for the 2-8 October period do not 
indicate any significant weather activity. A low pressure 
center was present southeast of the San Diego area early in 
the week, before moving farther southeast. Weather condi- 
tions were subsequently dominated by a strong offshore ridge 
which persisted throughout the week. 

Data were collected only on the 11th and 12th of 
October during the second week of TOWARD. Wind directions 
were from the northwest on the 11th and windspeeds were low. 
On the 12th, winds were also out of the northwest, but had 
much higher velocities. A maximum of about 9 m/ s was 
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(b) 



Figure 3.1 Surface analyses 1800 GMT 2 (a) 
and 3 (b) October 1984. 
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Figure 3.2 Surface analyses 1800 GMT 6 (a) 
and 7 (b) October 1984. 
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recorded on that date at 1550. Typical sea breeze patterns 
were also apparent in the data, particularly on the 12th. 
Weak easterly winds were present in the early morning. The 
winds shifted to the northwest near noon and increased 
rapidly to a maximum. These values began to drop off at 
approximately 1600 PDT. 

The synoptic charts from this period show a frontal 
passage through the San Diego area on the morning of 12 
October. This event cannot be correlated with the strong 
winds of the 12th, because the front was already several 
hundred kilometers to the southeast when the winds were 
recorded. Synoptic charts do not indicate the presence of a 
strong, large-scale pressure gradient associated with the 
trough. Thus the most likely cause of these high winds was 
strong diurnal heating and subsidence following the passage 
of the front. 

Data collection for the third week of the experiment 
began at 0900 PDT on October 15th. A definite wind shift can 
be seen on the 15th accompanied by sharp fluctuations in 
windspeed. Synoptic charts for that date indicated the pres- 
ence of a stationary front almost directly over the San 
Diego area at 1000 PDT. This front passed through the area 
at a rate of only 2 m/ s . Its slow rate of movement caused it 
to influence weather conditions in the San Diego area 
through the 16th. 

Wind shifts plus strong increases and decreases in 
windspeed were also apparent on the 16th and 17th. Synoptic 
charts for these dates show a front moving rapidly through 
the tower area at a rate of over 10 m/s. This front was well 
to the southeast within six hours . 

October 18th provided a more typical weather pattern 
with northwesterly winds and a strong increase in windspeed 
between 1000 and 1700 local due to the diurnal heating 
pattern. Synoptic charts show a strong high pressure ridge 
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Figure 3.3 Surface analyses 1800 GMT 11 (a) 
and 12 (b) October 1984. 
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(a) 




(b) 



Figure 3.4 Surface analyses 1800 GMT 15 (a) 
and 16 (b) October 1984. 
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(b) 



Figure 3.5 Surface analyses 1800 GMT 17 (a) 
and 18 (b) October 1984. 
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building to the southwest, with a weaker ridge to the 
northeast , 

Winds were westerly on the morning of the 20th with 
a shift to an easterly flow near midnight. Velocities were 
moderate and fairly constant, with values decreasing in late 
afternoon. 

The first observations for the week of 22-29 October 
began at 1100 on October 23rd. A strong increase in wind- 
speed can be seen, along with a northwesterly wind direc- 
tion. On the afternoon of the 25th, windspeeds again showed 
an increase. The wind direction was southerly in the late 
morning and shifted to the northwest by early afternoon. 
Wind speeds recorded on the 27th were moderate with little 
diurnal change. This pattern continued on the 28th, although 
the winds again originated from the northwest. 

The fifth week of the TOWARD experiment was from 29 
October to 5 November. Windspeeds were moderate on the 29th 
and 30th and were from the northwest. A wind shift and 
increase in windspeed began at approximately 0900 PDT on the 
31st. Evaluation of synoptic charts for 31 October indicated 
the presence of a large-scale pressure gradient associated 
with an oncoming front. On November 2nd, the winds were 
again out of the northwest and San Diego returned to the 
typical sea breeze pattern. 

The development of a strong high pressure center 
over the western United States produced the low windspeeds 
and shift in direction evident on 4 and 5 November. 
Windspeeds ranged from approximately 2 to 6 m/ s . 

Data collection for the final week of Phase I 
concluded on 7 November. Winds on the 5th were light and 
southwesterly. By the 7th, the typical sea breeze pattern 
again prevailed. Windspeeds were low and out of the west in 
early morning, but shifted to the northwest by early after- 
noon. A maximum value of about 5 m/s was recorded at 1105 
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(a) 




(b) 



Figure 3.6 Surface analyses 1800 GMT 20 (a) 
and 21 (b) October 1984. 
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(a) 




(b) 



Figure 3.7 Surface analyses 1800 GMT 23 (a) 
and 25 (b) October 1984. 
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(a) 




(b) 



Figure 3.8 Surface analyses 1800 GMT 27 (a) 
and 28 (b) October 1984. 
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Figure 3.9 Surface analyses 1800 GMT 29 (a) 
and 30 (b) October 1984. 
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(a) 




(b) 



Figure 3.10 Surface analyses 1800 GMT 31 (a) 
and 2 (b) November 1984. 
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(b) 



Figure 3.11 Surface analyses 1800 GMT 3 (a) 
and 4 (b) November 1984. 
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PDT followed by steadily declining values over the remainder 
of the sampling period. 
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(b) 



Figure 3.12 Surface analyses 1800 GMT 5 (a) 
and 7 (b) November 1984. 
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2. Phase II 



Data collection for Phase II of the TOWARD experi- 
ment began on 12 March 1985. It can be seen from Fig. 3.14 
that winds on the 12th and 13th were light to moderate and 
were primarily westerly. On March 15th and 16th winds were 
again moderate, but had a more typical northwesterly flow. 

Synoptic charts for this period do not indicate any 
significant weather activity. Conditions were dominated by 
an offshore high pressure center throughout the week. 

On March 18th, synoptic charts depict a tight pres- 
sure gradient resulting from the passage of a low pressure 
center through the San Diego area. The result was a sharp 
increase in windspeed. A maximum value of roughly 10 m/s was 
recorded at 1550 on that date. During the remainder of the 
week, the development of a strong offshore ridge reduced 
windspeeds and produced a northwesterly flow. 

The final week of the TOWARD experiment covered 25 
to 28 March. A gradual wind shift accompanied by light to 
moderate winds can be seen on the 25th. Synoptic charts for 
that date show a front passing through the San Diego area 
between 1000 and 1600 local time, at an approximate rate of 
4-5 m/s. As it moved westward, the low deepened and the 
frontal band remained southeast of San Diego through 1000 on 
the 26th. The winds reflect this frontal movement. A shift 
in direction to a more typical northwesterly flow occurs in 
late morning and a sea breeze pattern is apparent by late 
afternoon. The maximum wind speed recorded was about 7 m/ s 
at 1715. 

Winds on the 27th were out of the west - southwest and 
had generally higher magnitudes that those recorded the 
previous day. These values ranged from a mimimum of about 3 
m/s to a maximum of just over 8 m/ s at 1015. This maximum is 
supported by synoptic charts for 1000 that depict a strong, 
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(b) 



Figure 3.13 Surface analyses 1800 GMT 12 (a) 
and 13 (b) March 1985. 
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(a) 




Figure 3.14 Surface analyses 1800 GMT 15 (a) 
and 16 (b) March 1985. 





(a) 




Figure 3.15 Surface analyses 1800 GMT 18 (a) 
and 19 (b) March 1985. 
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(a) 




(b) 



Figure 3.16 Surface analyses 1800 GMT 20 (a) 
and 21 (b) March 1985. 
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(a) 




(b) 



Figure 3.17 Surface analyses 1800 GMT 22 (a) 
and 24 (b) March 1985. 
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large-scale pressure gradient associated with a deep low 
center over Vancouver, Canada. This system quickly moved 
southeast, where the low center deepened over Arizona and 
New Mexico for the next several days. The result was a tight 
pressure gradient which produced the high wind values 
collected on 28 March. On that date, a maximum wind value of 
over 12 m/ s was recorded at 1415. 
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(b) 



Figure 3.18 Surface analyses 1800 GMT 25 (a) 
and 26 (b) March 1985. 
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(a) 




(b) 



Figure 3.19 Surface analyses 1800 GMT 2.1 (a) 
and 28 (b) March 1985. 
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IV. DRAG COEFFICIENT RESULTS 



A principal motivation of this thesis is to examine the 
variation of the sea-surface drag coefficient with varying 
wind and wave states. Since windspeed magnitudes were gener- 
ally low to moderate (i.e. < 10-11 m/s) during both phases 
of TOWARD, the primary emphasis will be placed on Cjjj^ varia- 
tions as a function of wave state. 

The question of whether the neutral drag coefficient is 
constant or increases with windspeed, wave state and fetch 
has received considerable study in recent years. Garratt 
(1977) summarized sea-surface drag coefficients from 17 
experiments (based on eddy correlation and profile methods) 
published between 1967 and 1975. Observations of Cjjj^ (repre- 
sentative of a height of 10 m) plotted as a function of 
showed that Charnock's relation (Eqn. 26) was valid over a 
windspeed range of 4 to 21 m/ s . 

Large and Pond (1981) compared measurements of the 
momentum flux obtained by both Reynolds flux and dissipation 
methods. They reported that a comparison of both methods 
showed excellent agreement, on average, for windspeeds of 
4-20 m/s. Large and Pond also plotted neutral drag coeffi- 
cients against U^q- Although they found a slight increase in 
the average with windspeeds below 10 m/s, they felt that 

this relationship could be adequately described by a 
constant, where: 

lO^Cjjjj = 1.14 4 <U;l0'^^° 

Above 10 m/s, the values of demonstrated a more rapid, 

linear increase up to windspeeds of 26 m/ s . 

Smith (1980) used the eddy correlation method to obtain 
windspeed and heat flux measurements over the windspeed 
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range of 6-22 m/ s . He noted that the Charnock relation also 
held for neutral sea-surface drag coefficients measured for 
long fetch. These values were clearly seen to increase 
with increasing windspeed. Several points were not 
adequately explained by the Charnock hypothesis, however. 
These included: 

1. The dependence of the drag coefficient on windspeed 
was greater than predicted if the data at low windspeeds 
were not included. 

2. If the data in Garratt’s summary are viewed sepa- 
rately at windspeeds both above and below 10 m/s, the drag 
coefficient also appears to be nearly constant at windspeeds 
below 10 m/s. Smith's average Cjjjj value for low windspeeds 
was (1.11)10 . At higher windspeeds, this coeffcient tends 
to increase rapidly as a function of windspeed. 

Thus, all three studies indicate a relatively constant 
*^DN "va-lue at windspeeds below 10 m/s. It follows then that 
variations in the neutral drag coefficient at low to 
moderate windspeeds must be due to changes in the sea state. 

Garratt did not explore the behavior of with respect 
to fetch or sea state. He felt that current quantitative 
information was insufficient to adequately determine the 
dependence of the drag coefficient on factors other than 
windspeed. 

Large and Pond found the to be largely independent 
of fetch. For both long and short fetches, the mean drag 
coefficient remained at 1.14 x 10 for windspeeds less than 
10 m/s. 

Smith's comparison of with fetch showed higher drag 
coefficients at short fetches than those at longer fetches. 
This is to be expected since at short fetches actively 
growing waves absorb momentum from the wind while waves at 
longer fetches tend to be more in equilibrium with the wind 
field . 
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Smith also considered the variation of the neutral drag 
coefficient with sea state. A linear regression of vs. 
rms-wave height showed an increase in with an increase 
in wave height, although the correlation coefficient was 
rather low (.58). 

More recent studies have included a drag coefficient 
model discussed by Geernaert (1985a). Geernaert found the 
model exhibited considerable sensitivity to fetch. The 
influence of fetch on the model calculations was due to the 
wave spectral variations incurred by the dependence of the 
Phillips coefficient on the wave saturation parameter, 
Co/Uvf, where Co represents the phase speed. Geernaert 
observed that high drag coefficients often resulted from 
fetch limited conditions. 

Wu (1985) also looked at fetch as influenced by the 
Phillips spectral coefficient. He found larger spectral 
coefficients at shorter fetches than those at longer 
fetches . 

Measurements of the neutral drag coefficient over the 
sea by all methods often yield considerable scatter and the 
data must be carefully edited. Several types of data have 
been removed from this analysis based on dissipation 
methods. They include: 

1. Data collected during rain or fog conditions. Water 
droplets striking the hot film result in latent heat 
release. This sudden cooling induces a large jump in 
voltage, producing spikes in the data. In the inertial 
subrange particularly, noise from droplet contact contami- 
nates the high frequencies and results in very large drag 
coefficient estimates . 

2. Periods of light winds also produce erratic 
values as a result of convection off the tower mast . 

3. Flow distortion is a problem with light or moderate 
winds, particularly from the south. The best direction for 
data collection is from west to north. 
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4. Radio transmissions, particularly in the autumn, 
also posed problems. The hot films serve as antennae for 
transmitted energy, occasionally producing unreliable Cjjjj 
values . 

Analysis was based on plotting as a function of U^^q, 
as in the earlier studies by Garratt, Large and Pond, and 
Smith (Fig. 4.1). The data were then separated into two sets 
in order to compare short and long fetch conditions. Short 
fetch waves occurred when winds originated from north to 
east (Fig. 4.2), while long fetch conditions were approxi- 
mated during west to northwesterly winds (Fig. 4.3). 




TOWARD 



Figure 4.1 vs U^q for long and short fetches. 

*^DN again plotted against U-j^q for each fetch condi- 
tion. The mean and standard deviation for each 1 m/ s 
velocity band was calculated and can be found in Table III. 

Table III indicates that higher drag coefficient values 
did result from short fetch conditions. Mean values for 
both the 2-3 and 3-4 m/ s bands under short fetch conditions 
significantly exceeded those for long fetch. 
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Figure 4.2 Cjjjj vs U^q for short fetch conditions. 
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Figure 4.3 vs U^q for long fetch conditions. 

*^DN ■'^3.1ues for long fetch conditions also varied and 

as earlier predicted. Magnitudes were 



were not constant 
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TABLE IV 



^DN * ^10 Long and Short Fetches 

SHORT FETCH CONDITIONS: 360-115 DEGREES 



1 . 0-2 . 0 m/s 

(1 point) 

value=5 . 027 
s . d . = 

LONG FETCH 

1. 0- 2.0 m/s 

(1 point) 

value=2 .119 
s . d . = 

4. 0- 5.0 m/s 

(79 points) 

mean= .955 
s . d . = . 442 

7 . 0- 8 . 0 m/s 

(14 points) 

mean= .856 
s . d . = . 463 



2. 0-3.0 m/s 

(8 points) 

mean=2 . 017 
s.d.=.851 

CONDITIONS: 250-3 

2 . 0- 3 . 0 m/ s 

(53 points) 

mean=l . 485 
s.d.=.828 

5 . 0- 6 . 0 m/s 

(63 points) 

mean= .980 
s . d . = . 444 

8. 0- 9.0 m/s 

(8 points) 

mean= .939 
s . d . = . 254 



3. 0- 4.0 m/s 

(21 points) 

mean= 1.656 
s .d. = .407 

DEGREES 

3 . 0- 4 . 0 m/ s 

(81 points) 

mean= 1.432 
s . d . = . 850 

6 . 0- 7 . 0 m/ s 

(26 points) 

mean= .896 
s.d.=.377 

9.0- 10.0 m/s 

(4 points) 

mean= 1.413 
s .d. = . 104 



10.0-11.0 m/s 

(1 point) 

value= 1.479 
s . d . = 



fairly high for the 1-4 m/s bands, dropped off between 4-8 
m/ s , and subsequently increased between 8-11 m/s. 
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Overall, the dependence of the drag coefficient on wind- 
speed was found to be insignificant above 4 m/ s , based on 
scatter in the data as well as the large standard deviation 
at each windspeed interval. Below 4 m/ s , each there was an 
apparent increase in the Cjjjj with a decrease in windspeed. 
These values are somewhat suspect, however, due to the 
uncertainty in the quality of the data at very low 
velocities . 

Daily time series of wind stress were also examined, 
with both Cjjjj and U* plotted against time. Lack of data, low 
winds or strong stability prevented a close examination of 
the wind stress in some instances. Several case studies are 
described below, however. Appendix A contains figures rele- 
vant to the discussion. 

1. Frontal Cases: 15, 16-18 October and 26-29 March 

Fig. A. 6 shows the results when a front was passing 
the tower on 15 October. values ranged from .584 to 
2.058. Magnitudes can be seen to rise with rising winds 
until 0900, fall with decreasing winds between 0900 to 1030, 
and remain fairly steady with the winds near constant 
values. A sharp drop in windspeed at 1200 did not produce a 
corresponding drop in values. increased during this 
period of low winds to a maximum at 1500, and then dropped 
off sharply. Synoptically , Fig. A. 6 shows a wind shift 
between 1400-1500 for southerly to westerly winds with a 
corresponding increase in windspeed. values began rising 
shortly afterward and continued to increase up to the end of 
the sampling period depite a drop in windspeed beginning at 
approximately 1700. 

A similar frontal case is presented in Fig. A. 7 for 
16-18 October. Both U* and can be clearly seen to 
increase with rising winds to a maximum at 0600 on the 18th. 
A wind shift from southerly to northwesterly winds occurred 
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shortly after 6 am and several sharp drops in windspeed took 
place before the wind leveled off at approximately 6 m/s. 
The drag coefficient shows a similar trend but does not 
reflect the rapid decreases evident in the windspeed between 
0500 to 0800. 

On 26-29 March, passage of a front as well as south- 
erly winds produce high values with considerable scatter 
(Fig. A. 31). 

In each of these frontal cases, increased just 
before passage of the front. This trend is consistent with 
an earlier study by Denman and Miyake (1973) in which drag 
coefficients were seen to increase on the leading edge of a 
storm, and then decrease slightly or remain constant. They 
attributed this behavior to the nature of the wave field. 
Geernaert (1985b) also observed higher drag coefficients 
preceding an approaching front. He postulated that wave 
fields far ahead of and behind cold fronts are semi- 
independent of one another and are generally in equilibrim 
with local winds. Waves generated behind the front, often in 
a direction perpendicular to the baroclinic zone, can under 
certain circumstances travel through the front and interact 
with the wave field just ahead of it. Since the waves in the 
warm sector often run parallel to the local winds, the 
convergence of the two fields produces a much more energetic 
sea state just ahead of the front. The result would be 
higher estimates. 

2. Sea Breeze Cases: 18, 27, 29 October and 13, 15, 18 

March 

The sea breeze effect is apparent in much of the 
data from the TOWARD experiment. On each of the above dates 
(Figs. A. 8, A. 12, A. 14, A. 22, A. 23 and A. 25), the rises 
sharply and then decreases slowly to near equilibrium. 
Described in terms of surface roughness, drag coefficient 
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values are known to increase with growing wave state. As the 
wave field becomes increasingly saturated with momentum, 
absorption of energy slows, resulting in drag coefficients 
that approach a smaller and more constant value. 

3. Shift from Long to Short Fetch Conditions; 3-5 

November 

Data from 3-5 November depict a shift from long to 
short fetch conditions. Long fetch conditions are present 
until approximately 2300 on 3 November, as indicated by 

northwesterly winds. The drag coefficient estimates (Fig. 
A. 18) are correspondingly low. With the shift to northerly 
winds at 2300, a sharp increase can be seen in ^DN- Although 
these values decline over the next three hours , they again 
increase due to the variable wind direction produced by a 
series of wind shifts on the 4th. 

4. Variable Winds: 7 November 

High values as result of variable wind direc- 
tions can be seen between 0900 to 1200 on 7 November (Fig. 

A. 20). Magnitudes become more constant as winds shift to a 

steadier direction. 

5. Long Fetch and a Decrease with Time: 16,19, 22 

March 

These examples indicate long fetch conditions due to 
the northwesterly winds, as well as a general decrease in 
Cdn with time (Figs. A. 24, A. 26 and A. 29). Since long fetch 
waves are known to be more steady state, they tend to 
absorb less momentum than would short fetch waves. 

Consequently, the drag coefficient estimates have lower 
magnitudes and tend to approach a constant value more 
rapidly than those seen in earlier examples. 
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In summary, case studies suggest a dependence on 

sea state as illustrated by both surface roughness and the 
degree of equilibrium between U^q and the waves (i.e. the 
amount of momentum entering the wave field for growth) . 
Basic to this hypothesis on wind-wave coupling is that 
is generally independent of U^q for windspeeds between 4-9 
m/s , and that the large amount of scatter is associated with 
systematic trends in sea state, such as frontal and sea 
breeze conditions. This scatter should be reduced once wave 
state is incorporated as an additional parameter in 
formulations . 
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V. SUMMARY AND CONCLUSIONS 



Surface layer measurements covering a wide range of 
windpeeds have been analyzed for both phases of the TOWARD 
experiment. These measurements have been correlated with 
relevant mesoscale and synoptic features. The results were 
then interpreted relative to the 

Neutral drag coefficient values for both phases of 
TOWARD showed no significant windspeed dependence in the 
range between 4 and 9 m/s. The mean value was determined to 
be ( . 94 i 0.4)10^, which is somewhat smaller than results 
obtained by Smith (1980) and Large and Pond (1981). These 
differences appear to be insignificant, however, due to the 
scatter present in all three data sets. Cj^jj values are 
higher for windspeeds less than 4 m/sec, but these values 
are suspect as low windspeeds often produce unreliable drag 
coefficient estimates. 

Definite trends were observed in Cj^jj with fetch and sea 
state. Short fetch conditions produced generally higher Cjjj^ 
values than those for long fetch. Trends in sea state such 
as frontal and sea breeze conditions were reflected in drag 
coefficient estimates as well. As with Large and Pond, 
estimates were seen to increase sharply just prior to the 
passage of a front. During sea breeze conditions, the 
values increased with time. As windspeeds and the wave field 
reached equilibrium, however, the drag coefficients 
decreased toward a more constant value. 
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APPENDIX A 

TOWARD DAILY TIME SERIES 



This appendix contains daily time series for wind speed 
and direction, air and sea surface temperatures, the neutral 
drag coefficient and U*, for both phases of the TOWARD 
experiment. These figures provide a reference for the 
discussion provided in chapters three and four. 
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Figure A.l 2 October 1984. 
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Figure A, 2 3 October 1984. 
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Figure A. 3 6 October 1984. 
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Figure A. 4 7 October 1984. 



75 





11 OCT TOWARD 1984 13 OCT 



(O 

3 

5 

-j 

ly 

O 

<n 

Ui 

>y 

s 

C3 

Ui 

a 



21.0 



23.0 



20.0 



ir.o 



14.0 



AIRTSMP ... 
SEA TEMP 




A. 




Figure A. 5 11 to 13 October 1984. 
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Figure A. 6 15 October 1984. 
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Figure A. 7 16 to 18 October 1984. 
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Figure A. 8 18 October 1984. 
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Figure A. 9 20 to 22 October 1984. 
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Figure A. 10 23 October 1984. 
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Figure A. 11 25 October 1984, 
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Figure A. 12 27 October 1984. 
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Figure A. 13 28 October 1984. 
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Figure A. 14 29 October 1984. 



85 




30 OCT TOWARD 1984 31 OCT 



Figure A. 15 30 October 1984. 



86 



X 

a 

O 

09 

O 2.0 



LARGE i ?0 NO 



^ X 

xXx 



^ 




31 OCT 



TOWARD 1984 



01 MOV 





Figure A. 16 31 October 1984. 
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Figure A. 17 2 November 1984. 
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Figure A. 18 3 to 5 November 1984. 
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Figure A. L9 5 November 1984. 
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Figure A. 20 7 November 1984. 
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Figure A. 21 12 March 1985. 
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Figure A. 21 13 March 1985. 
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Figure A. 23 15 March 1985. 
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Figure A. 24 16 March 1985. 
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Figure A. 25 18 March 1985. 
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Figure A. 26 19 March 1985. 
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Figure A. 27 20 March 1985. 



98 




WIND ai^EED (M/aeC) OEQHEEa CEiaiUa u. EPaiLON lO^Co 







LARGE 1 POND ... 


- 




X 


■ 




X 






X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



XX . 

X X XX 



00 

21 MAR 



12 

TOWARD 1986 



00 

22 MAR 





Figure A. 28 21 March 1985. 
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Figure A. 29 22 March 1985. 
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Figure A. 30 24 March 1985. 
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Figure A. 32 26 to 29 March 1985. 
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